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We study the effect of a uniform magnetic field on tfie dynamics of axions. In particular, we 
show that the Peccei-Quinn symmetry is explicitly broken by the presence of an external magnetic 
field. This breaking is induced by the non-conservation of the magnetic helicity and generates an 
electromagnetic contribution to the axion mass. We compute the magnetic axion mass in one loop 
approximation, with no restriction on the intensity of the magnetic field, and including thermal 
. effects. 
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(D ■ I. INTRODUCTION 

c/^ : 

On , Axions are the pseudoscalar particles predicted by the Peccei and Quinn (PQ) mechanism a quite elegant 

• attempt to explain the smallness of the CP-violation induced by the QCD Lagrangian (the so-called strong CP 

I , problem; for a review see, e.g., Ref. 0)- They emerge as Pseudo-Goldstone Boson modes associated with the global 

j>. ■ axial PQ-symmetry [/(1)pq, spontaneously broken at the PQ-scale fa- This last is also known as the PQ- or axion- 

^\ , constant, and characterizes all the axion phenomenology. More specifically the axion mass is given by the relation 

On • 



a^ : 

o ^ 
\o 

o : 

Oh. 



X 



mqcD/eV ~ 6.2 x lO^GeV/^, (1) 

(here QCD means that this is induced by interaction with gluons) while its interactions with fermions are measured 
by Qi i^i/fa, where rrii represents the fermion mass (e.g. me,rnN, ... for electrons, nucleons, etc.). In addition, 
though chargeless, axions interact with photons through the electromagnetic anomaly: 



g^l Ca^^\gayaFF, (2) 



where F is the electromagnetic field, Ff^" = (1/2) e^^^-^Fap its dual, and = ctem£, / {"^i^ f a) ■ Here, acm is the 
electromagnetic fine structure constant and ^ an order one, model dependent constant. Presently, the PQ-constant's 
allowed range is rather narrow 0] : terrestrial experiments together with astrophysical and cosmological considerations 
have in fact excluded all the values of fa up to 10^ GeV, and above 10^^ GeV. ^ 

Because of the axion-photon interaction term |(2J), the phenomenology of axions can be influenced by the presence 
of an external electromagnetic field. Probably, the most discussed result of this well studied interaction is the axion- 
photon conversion in an external magnetic field (see, e.g., Ref. Q), which presently seems the best way to detect 
axions by terrestrial experiments 0] . 

Another interesting consequence of interaction (j^J was discussed in Ref. 10] where it was shown that axions in a 
strong magnetic field undergo a finite mass shift. However, in that computation of the mass shift the temperature 
effects were not taken into account. This makes the result in Ref. though correct, not directly applicable to the 
axion phenomenology. In fact, since the axion-photon coupling is very small, the induced mass can be relevant for the 
axion phenomenology only in the case of very intense magnetic fields, which normally exist only at high temperatures. 
In this paper, we reconsider this effect accounting also for the thermal contribution. As we will show, the result given 
in Ref. [lO| is a very good approximation only for temperatures well below the electron mass T <C m^. For higher 
temperatures and sufficiently intense magnetic fields the thermal corrections are relevant. 



^ In the case of the hadronic axion [ij, a small window around fa ~ lO^'GeV can be also permitted (see also Ref. 0). 
^ Attempts to relax the above bounds on the PQ-constant arc discussed in Ref. whereas possibilities of different cosmological bounds 
are discussed in Ref. 0- 



2 



Interestingly, the resulting mass shift induced by the magnetic field is independent of the gluon dynamics and, 
therefore, does not vanish if the (QCD induced) axion mass 1^ is set to zero. In other words, the external electro- 
magnetic field itself is responsible for the breaking of the PQ-symmetry, playing a role similar to the gluon fields. We 
will discuss this point in Section 2 of this paper. We will show that the magnetic contribution is indeed due to a 
topological effect (just as for gluons) which is expected when the axion dynamics take place in an external magnetic 
field. Ultimately, we will associate this effect with the generation of magnetic helicity, whose important role in axion 
cosmology has already been pointed out (see Ref. ^3)- This property makes the magnetic mass shift a true mass term. 
Therefore, the terminology "mass shift" is perhaps inappropriate. In the following, we will refer to this contribution 
as the magnetic induced axion mass niB- It contributes to the total axion mass as 

m'^ ^m'^QCD+m%, (3) 

where toqcd is given in 

Because of this property, our result applies not only to axions but to any pseudoscalar particle which interacts 
with photons as in Eq. independent of its interaction with gluons. This is the case with arions 12], (almost) 
massless pseudoscalar particles which do not interact with gluons at all. On the other hand, concerning the axion 
phenomenology, we expect the contribution ttib to be especially relevant at temperatures above ~ IGeV, where the 
gluon contribution to the axion mass is strongly suppressed 15], or in regions of space and time where a very strong 
magnetic field is present. A more detailed analysis of the phenomenological implications of the magnetic induced 
axion mass is in progress. 

The paper is organized as follows: In Section 2 we comment on the origin of the magnetic induced mass term. We 
then compute the magnetic induced axion mass in Section 3, discussing the zero temperature limit as well as the finite 
temperature result. Finally, in Section 4 we summarize and give our conclusions. 



II. PECCEI-QUINN SYMMETRY BREAKING BY AN EXTERNAL MAGNETIC FIELD 

In this section we want to discuss some features about the origin of the axion mass generation due to the external 
magnetic field. This discussion is semi-qualitative and based only on symmetry considerations. The actual compu- 
tation of the induced electromagnetic mass term is devoted to the next section. However, here it will be shown the 
relevance of the magnetic helicity for axion dynamics (see also Ref. [HI 

In the PQ-mechanism the axion is defined as the phase of a field (or, more generally, a linear combination of the 
phases of different fields). The Lagrangian is constructed so that it is invariant under the PQ-symmetry [/(1)pq. At 
pure classical level the PQ-symmetry is exact and can be read as the invariance of the effective action Ses[a{x)] [see 
Eq. for a constant shift of the axion field a(x) a{x) + c {c = c' fa, with c' a constant). This symmetry requires 
the axion interactions to be of the form J^^ (where J^^ is a conserved current), and prevents the generation of the 
mass term rn^a^ . If quantum mechanical effects are taken into account, it results d^J^ ^ 0, and the interaction terms 
with the gluons [a/ fa) GG and the electromagnetic field {a/ fa) FF emerge. Both GG and FF are total derivatives. 
However, because of the non-trivial topology of the gluon fields, the first one is responsible for the axion mass term. 

Following Shifman, Vainshtein, and Zakharov ,4] we can state the problem better defining the quantum effective 
action for the axion field, obtained by integrating away all the fields but the axion: 

g-Scft[a(a:)] = j q-S[o.{x),4>] ^ 

where is a collective name for the gluon, photon and fermion fields. Since the derivative terms like are left 
invariant for a + c, it follows that: 

g-Srft[a(a:)+c] ^Jj^^^-Sla{x),<t>] + cgJd-^xGG + cg^^Jd^xFF ^g-j 

{g is the coupling with gluons). If we neglect the external magnetic field we can also neglect the term J d'^x FF since 
we know that the topology of the electromagnetic field is trivial. However, the / d'^x GG term cannot be neglected and 



^ Since arions do not interact with gluons, they are not affected by the color anomaly. Therefore, they do not have a QCD-induced mass 
term as in Eq. The phenomenology of arions in an external magneti c fie ld has been re con sidered recently as a tentative to analyze 
the problem of supernova dimming in terms of photon-arion oscillations |13| (see also Ref. 
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is therefore responsible for the breaking of the PQ-symmetry, Scs[a{x)] — > Scs[a{x) + c], and the origin of the axion 
mass. Let us now suppose that an external magnetic field is present in addition to the radiation one, and reconsider 
the / d'^xFF term. We use the relation FF = ^^'K^, where = (e^„/37^"5^^^) = {I<o, -K) and Ai" = (Aq, A) is 
the total electromagnetic field, which includes the external field. Neglecting the contribution of A at the surface of 
the spatial integration volume V, which implies Jy (PxdrvJi. = 0, we get 



*2 

d^xF^.F^" ^HB{t2)-HB{ti), (6) 



where 



HB{t) = ^ I d^xA- V X A, (7) 



is known as magnetic helicity or Abelian Chern-Simons term (see, e.g., Ref. ^3)- The conclusion is therefore that 
the term J d'^x FF cannot a priori be neglected if the helicity of the system is not conserved. 

Interestingly, it was shown a few years ago \\^ that the dynamics of the axion field (indeed of any pseudo- 
scalar field) in an external, uniform magnetic field necessarily induces a change of the magnetic helicity. Therefore, 
the quantity HB(t) is not constant and, consequently, J d'^xFF ^ 0. The non- vanishing of this term implies the 
non-invariance of the classical action for a constant shift of the axion field, S[a{x) + c] — S[a{x)] = cga-y J d'^xFF cx 
AHb 7^ 0. This translates in a mass term for the axion field, if also the quantum effective action is non-invariant 
Scs[a{x) + c] S'efi[a(x)]. This is the case if all the charged fermions are massive. On the other hand, if at least one 
charged fermion ip (with electric charge g^) is massless, we could redefine the fermion field as 

^P^^' = e"^''" V', = -0 e'"^" , (8) 

with a constant. This would introduce a change in the measure 

dlpdip ^ dWd^p' ^ e^°''''l^'^^''^^ dlpd^, (9) 

(a' is a constant proportional to a) and no changes in the fermion mass matrix. Therefore, it would always be possible 
to choose a so that (c + a' ) J d'^x FF = 0. In other words, if one charged fermion is massless it would be possible to 
reabsorb J d^x FF by an unphysical chiral rotation of that fermion, and then the axion would not get a mass (this is 
analogous to what happens to the standard QCD-induced mass if a quark is massless). 



III. COMPUTATION OF THE MASS 



In this section we compute the effective, magnetic-induced, axion mass. Technically, the problem consists in the 
computation of the vacuum polarization effects produced by the external (classical) electromagnetic field. Diagram- 
matically we consider the temperature dependent vacuum polarization tensor 11^^ with one loop quantum corrections 
and with an arbitrary number of interactions with the external field. 

For what concerns our conventions we use natural units with h = c = \ and the Lorentz-Heaviside convention so 
that = e^/(47r) ~ 1/137. The metric tensor is chosen as — diag(l, —1, —1, —1) and we adopt the often used 
notation 

fc|| = (fco, 0,0,^3), 

fc± = (0,A:i,fc2,0), (10) 

for the axion momentum = (fco,k). In addition, we assume the external magnetic field Bext, with intensity 
|Bcxt| — B, to be directed along 2, Bext = Bz. The total electromagnetic tensor is F^^'^ = F^^^ + F'^'^ , with 
pfiv _ g^j^^ — di,Afj_ accounting for the quantum electromagnetic field Afj_, and i^^xt for the classical external one. 

In order to compute the magnetic induced axion mass we need to consider the exact axion propagator in the external 
field, A(A:), which is related to the a ^ a transition amplitude in the external field, M{k), by the equation A(fc) = 
Ap{k) + AF{k)M{k)A{k), where Apik) is the free axion propagator. The above equation is shown diagrammatically 
in Fig. 1. Its formal solution is A"^ = k^ + rnqQ^y — M (k) . The magnetic axion mass can be calculated in the standard 
way taking the limit k ^ and solving the equation = —'ReM{mB)- 
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FIG. 1: Diagrammatic representation of the equation A{k) — Apik) + AF{k)M{k)A{k) for the exact axion propagator in the 
external magnetic field Boxt- The sources of this are denoted by crosses. The dashed lines represent the free axion propagator, 
while the double dashed ones represent the exact propagator. The internal waved line stands for the exact photon propagator 
in the external magnetic field. 



We are considering here only the interaction of the axion with the photon. Of course, also the axion-fermions 
interactions contribute to the axion propagator, but the resulting mass shift is small and disappears in the limit of 
™QCD Therefore, in the following, it is assumed that M{k) is induced only by photons (as in Fig. 1). 

Formally, the exact photon propagator is given by the Dyson relation 

V~\k) ^ D^^k)~Il{k), (11) 

where D^^^, — {^\D{k)\i/) — k^g^v — Qkfiki, is the free propagator and C, is the gauge parameter. The polarization tensor 
n is diagonalizable and its 4 eigenvectors le*^'-*), i — 1, .., 4 satisfy the orthonormality (e^'^le^^-*) — Sij and completeness 
|e^'-*)(e'^*^| = 1 relations Because of the transversality condition, A: • 11 = 11 • = 0, the polarization tensor has 
a zero eigenmode directed along k. To keep the same notation as in Ref. we label this as the fourth eigenvector 
e^^'> = {fi\e^^^) = kf,/\k\. So we find 

n = j2>^^\e^'^){e^'^l (12) 



where we have indicated with 



A, = (e«|n|e«) = ef)n^,e(^) (13) 



the eigenvalue associated with |e^*^). Of course, it results A4 = 0. Using the orthonormality and completeness relations, 
it is easy to invert Eq. Hll(l to get a formal expression for the exact photon propagator 

i—l 

In an external magnetic field the relevant term of the axion-photon interaction Q can be conveniently written as 

^a-f — ga'jO- Fl^^^ d^A^. (15) 

Observing that the vector e^^-* k'^ F^;^^ / {\k\\\B) directed along Bcxt is an eigenvector because of the azimuthal 
symmetry induced by Boxt j this interaction takes (in momentum space) the particularly simple form 



Caj = -igaj\k\\\BaA^e';^\ (16) 
Using the orthogonality between the eigenvectors together with Eq. (|14() . we find 

where the eigenvalue A3 is given by Eq. l(T^ . In the hmit k — > 0, it results A3 = 1133. Let us define the function V{ko) 
such that 



n3#5^k^^ (18) 

Stt kf 
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where we have introduced the factor acm/37r for convenience. It results 

M(fco)^- , , ^^^'T.-^,,. - (19) 
1 + (Q;om/37r)>'(A;o) 

The one loop polarization tensor at finite temperature in an arbitrary magnetic field was computed using the imaginary 
time formalism in Ref . |2fij | . In the imaginary time formalism, the (Euclidean) photon energy is quantized and expressed 
in terms of the so-called Matsubara frequencies: fc4 = —ik^ ujn = 27rnT, with n integer. For simplicity, we consider 
only the contribution of the electrons in the fermion loop since this is the leading one as electrons are the lightest 
charged fermions. In any case, the contribution from the other fermions can be easily included. The electron energy 
in the fermion loop is quantized as well, p4 = — ipo = (2/ + l)7rr, with I integer. The function V{lu„) is given in 

Ref. 20]. However, for our purpose, the computation of the effective axion mass can be simplified observing that in 
all the phenomenologically relevant cases the axion mass is always smaller than the fermion masses ms <C We • This 
assumption will be verified a posteriori. Therefore, we can safely consider the limiting value V = V^LOn/m,, — > 0) in 
place of V{uJn)- So finally for the magnetic- induced axion mass we find 



^ 1 + (a„n/3^)P 

From Ref. |20|| we have 



/— — oc 





u coth u — 1 j 







(20) 



(21) 



where u = eBu. Observe that the above expression is convergent and the last term in the right hand side of Eq. H21() 
represents the contact term which removes the ultraviolet divergence (u 0). The normalization chosen gives uib — 
in the zero temperature and zero magnetic field limit. 



Ilia. Low Temperature Limit 

Before a complete analysis of expression (I21|) let us consider its low temperature limit Vq. Actually, the computation 
of Vq is straightforward since in the limit T ^ the sum over the Matsubara frequencies of the electrons can be 
substituted by an integral over the (Euclidean) electron energy: Tj^i ^ (27r)^^/dp4. The result, 

f°° du 2 

= lim P = / — e-™= [ucothw - 1] , (22) 

is in agreement with the existent literature (see Ref. |lOl|\ To have a better comparison with the non-zero temperature 
result it is convenient to introduce the dimensionless variable 

B eB 
K 

so that 



(23) 



/•oo 

Vo= dse-"'^ 
Jo 



coth s 

s 



(24) 



(The parameter Be = m^/e ~ 4.4 x IQ-'^^G is known as the critical or the "Schwinger value" for the magnetic field.) 
A good approximation of the above integral (see Fig. 2) is given by the following expression: 

7^0 ^ A coth A- 1. (25) 

In particular, the limit A ^ 1 (S <C Be) gives Vq ~ A^/2, whereas in the opposite limit X ^ 1 {B ^ Be) wc find 
T'o — A. Taking into account Eq. H20I) . this last limit gives 



2 dafB 

l-t-(aem/37r)(B/S,) ' 



(26) 

± -r yucm/ on jy-u/ ue) 
which for B/Be^ 37r/acm ^ 10^ gives in turn the result 



presented in Ref. and used in Ref. [2l| . 
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FIG. 2: The continuous line represents the function Po in Eq. (|24j, while the dotted line is its approximate profile Eq. ||2 



Illb. Finite Temperature Results 

In order to include the thermal effects it is necessary to calculate the sum over the Matsubara frequencies in Eq. H21|l : 



E 



e---. ( 1 _ _ u^tu 



This sum can be conveniently expressed in terms of the the Jacobi theta function of second kind [22 

oo 

d2{w, q) = 2q^/^ cos[(2n + l)^. 



n=0 



Indeed, observing that 



/ — — OO 



E 

^ — — oc 



du 



E 



we have 



Ee 

/— — oo 



To proceed further, we define the function 



2 d 

3 du J 



(28) 



(29) 

(30) 
(31) 

(32) 



(33) 



whose graph is shown in Fig. 3. Taking into account Eqs. I|32() and (|33|l . we re-cast Eq. (|21|1 in the following form: 



where we have defined 



V= dse-'/^ 
Jo 



t = 



J^(st^)coths- - 
s 



ttT 



(34) 



(35) 
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FIG. 3: The function T{x) defined in Eq. CT . 



Equation 134|l is very similar to Eq. H24|l. except for the presence of the function T{st^). It is important to notice 
that the temperature enters only through !F{sP), and only as sP. Let us comment this result. First of all we see 
from Fig. 3 that J-'{x) ~ 1 for x sufficiently smaller than 1. In addition, in Eq. (|34|l the parameter A has the role of a 
cutoff. Therefore, in the relevant integration region, < s < few x A, it results st^ < few x Xt^ = {TiT/mef'. Then, 
if T ^ TOe, we can set ~ 1 and consequently we find the result (|24|l . On the other hand, for T ^ and B < T^, 
we find the asymptotic expansion of V: 

V ~ 7 + lnc-lnf —\ , (36) 

where 7 ~ 0.577 is the Euler's constant, c ~ 2/3 is a numerical constant. If i? > T^, the above equation is still valid 
with c a slowly decreasing function of BjT^ (see Fig. 4). 

Our numerical results are shown in Fig. 5 which represents V given in Eq. (|34|l as a function of the normalized 
temperature T/rrie, for different values of the parameter A — B/Bc- Also represented is the asymptotic expansion 
Eq. H36() . A look at Fig. 5 shows that for a large range of values for the external magnetic field and temperatures, 
the term (acm/37r)7' is much smaller than 1. In particular, this is the case for high temperatures (T 3> me) whatever 
the intensity of the magnetic field is, or for low temperatures (T ^ me) and magnetic fields below the critical value 
[B < Be). In these cases, from Eq. ()2()|l we get that the axion mass is well approximated by 

ms ^ gajB. (37) 

In order to quantify the thermal effects on the axion mass, we compare the above result with the axion mass as 
calculated in the zero-tcmperature limit and for large magnetic field (see Section Ilia). From Eqs. (|27() and H37() we 
get 

mliT,B) B_ 

m|(0,S) ~ 3^ Be ■ ^ ' 

Therefore, for T 3> me and B ^ Be, we find that the axion mass is strongly influenced by temperature effects. In 
particular, thermal effects give rise to in increment of the axion mass with respect to the zero temperature limit. 

Before closing this section, a comment is in order. The complete result in Eq. 13()|l is somewhat puzzling, since it 
gives a singular limit when me — > 0, whereas we would expect m^ ^ in the chiral limit, as discussed in Section II. 
However, it is easy to understand that this limit can not be considered in our discussion. In fact, the approach that 
we have considered is based on the perturbative expansion of Fig. 1, and loses clearly its validity if M{k) is arbitrarily 
large, as in the case of very small electron mass. Indeed, our result is based on an even stronger assumption, which 
is that the electron mass is always much larger than the axion mass. Therefore, the result (|2(JII is valid only if the 
inequality 



(39) 
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FIG. 4: The function V defined in Eq. (gU versus the temperature for difi'erent values of the "normalized" magnetic field 
h = B/T'^. Dotted line: b = 100; Dashed line: h = 50; Long-dashed line: b = 10; Solid line: h = 1. For & < 5, all curves collapse 
into the fine = 1, which is well approximated by Eq. 13611 . 
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FIG. 5: Numerical representation of the function V given in Eq. 1)34^ as a function of the normalized temperature T/nie for 
difi'erent values of the parameter A = B /Be- Left panel. From the top to the bottom: A = 20, 15, 10, 1. Right panel. From the 
top to the bottom: A = 0.8, 0.6, 0.4, 0.2. The tick solid line in both panels represents the asymptotic expansion l|86^ . 



is verified, which is not the case in the chiral limit. However, for the standard value of the electron mass, the above 
inequahty is essentially always verified in any physical interesting case. In fact, even in the "worst" case P ^ 1 (which 
refers to the case T <C We and B Be, or to the case T ^ nie), the above condition translates in 



which is satisfied, even in the most conservative case fa = lO^GeV, by any cosmological and astrophysical magnetic 
field present in the universe below the electroweak phase transition ,23|| . 



IV. CONCLUSIONS 



The phenomenology of the axion in an external magnetic field has been deeply studied in the last few decades. 
Nowadays the relevance of an in depth understanding of the effects of axion interaction in an external magnetic field 
to astrophysics and cosmology is clear. In this paper we have further investigated this problem computing the axion 
mass induced by an external magnetic field. The origin of this phenomenon is related to the generation of magnetic 
helicity, induced by the axion itself. We computed the resulting induced mass in one loop approximation, accounting 
for the thermal corrections, and with no restriction on the intensity of the magnetic field. Our result is given in 
Eqs. (|20|l and H34|) (see also Figs. 4 and 5), and indicates that thermal effects are not relevant for temperatures 



9 



significantly below the electron mass. For T <C our result is indeed in agreement with the existent literature. 
However, for higher temperatures, the axion mass is strongly influenced by temperature effects if the external magnetic 
field is much stronger than the Schwinger value, B ^ Be ^ lO^'^G. In this case, thermal effects give rise to an 
increment of the axion mass with respect to the zero temperature limit, which is proportional to the magnetic field: 
ml{T,B)/mliO,B) ~ {aem/Sir) (B/B,). 

It is clear that the value of the magnetic field that could induce a sizable change of the axion dispersion relation is 
very high, and certainly well above the possibility of the present terrestrial experiments. However, in an astrophysical 
or cosmological contest the magnetic induced mass could play an important role. This topic will be the object of 
future investigations. 

We would like to thank P. Cea, V. Laporta, R. D. Peccei, G. G. Raffelt, and M. Ruggieri for helpful discussions. 
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